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ABSTRACT 

We present an overview of four ah initio axisymmetric core-collapse supernova simulati ons employing de- 
tailed spectral neutrino transport computed with our Chimera code and initiated from IWoosley & Hegeri 
(^007) progenitors of mass 12, 15, 20, and 25 Mq. All four models exhibit shock revival over ^ 200 ms 
(leading to the possibility of explosion), driven by neutrino energy deposition. Hydrodynamic instabilities that 
impart substantial asymmetries to the shock aid these revivals, with convection appearing first in the 12 Mq 
model and the standing accretion shock instability (SASI) appearing first in the 25 Mq model. Three of 
the models have developed pronounced prolate morphologies (the 20 Mq model has remained approximately 
spherical). By 500 ms after bounce the mean shock radii in all four models exceed 3,000 km and the diagnostic 
explosion energies are 0.33, 0.66, 0.65, and 0.70 Bethe (B = 10^^ ergs) for the 12, 15, 20, and 25 Mq models, 
respectively, and are increasing. The three least massive of our models are already sufficiently energetic to 
completely unbind the envelopes of their progenitors (i.e., to explode), as evidenced by our best estimate of 
their explosion energies, which first become positive at 320, 380, and 440 ms after bounce. By 850 ms the 
12 Mq diagnostic explosion energy has saturated at 0.38 B, and our estimate for the final kinetic energy of the 
ejecta is ~ 0.3 B, which is comparable to observations for lower-mass progenitors. 
Subject headings: neutrinos — radiative transfer — supernovae: general 



1. INTRODUCTION 



iColgate & Whita ([T966h originated the idea that core- 
collapse supernovae (CCSNe) may be neutrino driven. Be- 
cause realistically modeling CCSNe is a computationally 
complex task of tying together hydrodynamics, neutrino 
transport, general relativity (GR), nuclear burning, and a re- 
alistic equation of state (EoS), and since the many feed- 
backs and the nonlinearity of the process preclude an ana- 
lytical treatment, progress in the field has been closely tied 
to adv ances in computational resources and code develop - 
ment (iMezzacappal I2QQ5I: iKotake et al.l I2QQ6I: Uankal IMl*). 
Two decades ago advances of this sort made possible the first 
generation of CCSNe simulations in two dimensions (2D) 
Jllerantetal. 1994; Burrows et al. 1995; Janka & Miille^ 
|1996; Mezzacappa et al. 1998; Eryer & W arrenll200 2). These 
simulations, by allowing, among other things, neutrino- 
driven convection below the shock and continued accre- 
tion and shock expansion to coexist, gave explosions. 
In contrast, continually refined spherically- symmetric CC- 
SNe simulations failed to explode except f or the lightest 
progenitors (.Rampp & Janka.2000: .Mezza cappa 
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iLiebendorfer etaH I2QQTI I2QQ1 iBuras elaH l2QQ3h . Eurther 
investigation has shown the crucial importance of multi-D 
effects, including the Stand ing Accretion Shock Instability 
(SASI) (iBlondin et al.ll2003h . in shaping explosions, and in 
accounting for the many observed CCSNe explosion asym- 
metries (see, e.g., Janka 2012). 

While the first-generation 2D CCSNe simulations were 
limited by their use of gray neutrino transport or precom- 
puted spectral neutrino transport, current computing resources 
have permitted second-generation 2D CCSNe simulations 
with self-consistent, spectral neutrino transport (Buras et al. 
2003, 2006b a; Burrows et al. 2006, 2007; Bruenn et al. 2006; 
Bruenn et al. 20 09: Marek & Janka 2009; Suwaetal. 2010; 
iTakiwaki et al.l l2012: Miiller et al. 2012b, a), essential for ac- 
curately computing the neutrino-matter coupling at the heart 
of the neutrino-driven CCSN mechanism. Realistic multi-D 
ah initio (neutrino luminosities and spectra computed from 
detailed spectral neutrino transport rather than prescribed, pa- 
rameterized, or incomplete spectral transport) CCSNe simu- 
lations are essential for elucidating the details of the neutrino- 
driven supernova mechanism and fully accounting for super- 
nova observables, including explosion energies, neutrino and 
gravitational wave signatures, nucleosynthesis, compositional 
mixing, and neutron star recoils. The Garching group has 
published such simulations for the 2D case using the Ver- 
tex code, incorporating a complete set of neutrino opaci- 
ties, velocity dependent transport terms, and neutrino-energy- 
bin coupling in the neutrino transport, GR corrections to 
Newtonian self-gravity or approximate GR self-gravit y, and 
standard Eo Ss (Buras et al. 2006b, a; Marek & Janka| |2009|; 
iMiiller et al.l I2012bll3) . Their successful models either pro- 
duce weak explosions or were reported before their explosion 
energies saturated. 

The limited number of published ab initio CCSNe simu- 
lations with detailed neutrino transport, the quantitative dif- 
ferences we obtain relative to the Garching group, and the 



2 Bmenn et al. 



Table 1 

Model Summary Table 



Models 



Property 


rSiZ-WrlU/ 


ri ID- WrlU / 


oZU- WrlU / 


13 o< wixjnn 
oZD- WrlU / 


Progenitor mass, Mzams [Mq] 


12 


15 


20 


25 


Progenitor compactness, ^1.75 


0.263 


0.332 


0.461 


0.465 


Initial simulation domain radius [km] 


30000 


20000 


20000 


20000 


Simulation enclosed mass [Mq] 


2.08 


2.75 


3.51 


4.57 


Radial zone count 


512 


512 


512 


512 


Angular zone count 


256 


256 


256 


256 


Collapse time (initiation to bounce) [ms] 


266 


335 


462 


470 


Bounce, central density [lO-*^^ g cm~^] 


3.13 


3.45 


3.42 


3.29 


Bounce, central Ye 


0.240 


0.249 


0.248 


0.246 


Bounce, shock position [Mq] 


0.45 


0.46 


0.47 


0.47 


Time to growth in gain region mass, tgain [ms] 


256 


221 


195 


192 


Time to positive diagnostic energy > 0.01 B [ms] 


219 


207 


192 


197 


Time to mean shock radius of 500 km, ^500 [ms] 


236 


233 


208 


212 


Post-bounce epoch, tpb [ms] 


850 500 


500 


500 


500 


Post-bounce diagnostic energy, [B] 


0.38 0.33 


0.66 


0.65 


0.70 


Post-bounce mean shock radius [km] 


7245 3429 


3293 


3632 


3783 


Post-bounce shock deformation, dghock 


0.42 0.31 


0.52 


0.027 


0.31 


Post-bounce PNS rest mass [Mq] 


1.48 1.48 


1.66 


1.80 


1.88 



specifics of our findings (explosions more broadly consistent 
with observations) motivate us to present in this Letter an out- 
line of four simulations being performed with our Chimera 
code. These simulations span progen itors of mass 12-25 Mq 
with detail comparable to Marek & Jankal (|2QQ9|) . Here we 
describe some of the general features already manifest in our 
simulations and how they support the viability of the neutrino- 
driven CCSN mechanism. 

2. NUMERICAL METHODS AND INPUTS 

Chimera is a parallel, multi-physics code built specif- 
ically for multidimensional simulation of CCSNe. It is a 
combination of separate codes for hydrodynamics and grav- 
ity; neutrino transport and opacities; and nuclear EoS and 
reaction network, coupled by a layer that oversees data 
management, parallelism, I/O, and control. The hydrody- 
namics are evolved via a dimensionally- split, Lagrangian- 
Remap (PPMLR) scheme (Colella & Woodward 1984) as 
implemented in VHl (Hawleyetal. 2012). Self-gravity is 
computed by multipole expansion ( Miiller & Steinmetz 1995) 
replacing the New tonian monopole with a GR monopole 
(iMarek et al .1120061 Case A). Neutrino transport i s computed 
in the "ray-by-ray-plus" (RbR+) approximation (iBuras et al.l 
l2003h . where an independent, spherically symmetric transport 
solve is computed for each angular "ray" (^-zone). Neutri- 
nos are advected laterally (in the ^-direction) with the fluid 
and contribute to the lateral pressure gradient where p > 
10^^ g cm~^. The transport solver is an improved and up- 
dated version of the multi-group flux-limited diffusion trans- 
port s olver of Bruenn (1985), enhanced for GR (Bruenn et al. 
|2001|) . with an additional geometric flux limiter to prevent the 
over-rapid transition to free streaming of the standard flux- 
limiter. All 0(v/c) observer corrections have been included. 
We solve for all three flavors of (anti-)neutrinos with four 

coupled species: Ue, i^e, ^fir = {l^^.l^r}. I^fir = {z^/i,^r}, 

with 20 energy groups each for ae = 4-250 MeV, where 
a is the lapse function and e the comoving-frame group- 
center energy. The neutrino-matter interactions include emis- 
sion, absorption, and non-i soenergetic scattering on free nu- 
cleons (Reddv et af] Il998h with weak magnetism correc- 
tions (Horowitz 2002^; emission/absorption (electron cap- 
ture) on nuclei (Langanke et al.„.2003,) : isoenergetic scatter- 



ing on nuclei, including ion-ion correlations; non-isoenergetic 
scattering on electrons and positrons; and pair emission 
from e^e~ -anni hilation (Bruenn 1985 ) and n ucleon-nucleon 
bremsstrahlung (iHannestad & Raffeltl Il998h . Unlike the 
Garching group, we do not include the effective mass cor- 
rections at high density nor the i^e^e ^/ir^/ir process. 
We utilize the K = 220 MeV incompressibility version of 
the Lattime r & Swestyl (1199 Ih EoS for p > lO^^gcm"^ 
and a modified version of the ICoopersteinI (11985) EoS for 
p < 10^^ g cm~^ where nuclear statistical equilibrium (NSE) 
applies. To aid the transition to the XNet (iHix & Thielemann! 
1999) 14-species a-network (a, ^^C-^^Zn) used for the non- 
NSE regions, we have constructed a 17-species NSE solver to 
be used in place of the Cooperstein EoS for Ye > 0.46. An 
extended version of the Cooperstein electron-photon EoS is 
used throughout. 

During evolution the radial zones are gradually and auto- 
matically repositioned during the remap step to track changes 
in the radial structure. To minimize restrictions on the time 
step from the Courant limit, we "freeze" the lateral hydro- 
dynamics for the inner 8 zones during collapse, and, after 
prompt convection fades, we expand the laterally frozen re- 
gion to the inner 8-10 km. In the "frozen" region we set 
ve = and skip the lateral hydrodynamic sweep. The full 
radial hydrodynamics and neutrino transport are always com- 
puted to the center of the simulation for all angular rays. A 
more extensive description of Chimera is under preparation. 

Four non-rotating axisymmetric models (designated B12- 
WH07, B15-WH07, B20-WH07, B25-WH07, corresponding 
to progenitor masses of 12, 15, 20, and 25 Mq), and identical 
ID models, are initialized without artificial perturbation from 
the pre- supernova progenitors of Woosley & Heger ( 2007). A 
grid of 512 non-equally spaced radial zones covers from the 
stellar center into the oxygen-rich layers. The 2D models em- 
ploy 256 uniform- sized angular zones from to tt. 

3. RESULTS 

All 2D models experience similar evolution immediately 
after bounce. At a post-bounce time, tpb, of a few ms the 
shock breaks through the -sphere, initiating a rapid delep- 
tonization front whose inner edge advects inward with the 
fluid. An unstable lepton gradient results, which is stabilized 
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Figure 1. Evolution of the entropy (upper half) and radial velocity (lower half) for B12-WH07, with snapshots at fpb = 12, 90, 150, 200, 300, 400, 600, and 
800 ms. The scale grows in time to capture the expansion of the supernova Shockwave, but the colormaps remain constant. The radial velocity portion is omitted 
for the first two snapshots. (An animated version of this plot is available at http : / /astro . phys . utk . edu/activities%3achimera%3aseriesb) 
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Figure 2. Upper panel: Mean, minimum, and maximum shock radii versus 
tpb for all simulations. Lower panel: Mean shock radii for all simulations 
and their ID equivalents. 

by the positive entropy gradient laid down by the shock as 
it initially strengthens. The initial oscillations of the proto- 
neutron star (PNS) periodically strengthen and weaken the 
shock, causing entropy-unstable pockets and transient con- 
vective episodes behind the shock. By tpb ~ 12 ms an ex- 
tended unstable entropy gradient has developed behind the 
shock as it weakens, extending down to ~ 60 km and driv- 
ing a brief episode of convection (Figure [H 12 ms panel). 
This episode of "prompt" convection has little appreciable ef- 
fect on the shock, and by tpb ~ 40 ms convective activity has 
ceased. 

Up to tpb ~ 100 ms, the shock trajectories of correspond- 
ing ID and 2D models track each other closely (Figure O 
lower panel). Remarkably, the shock trajectories for differ- 
ent progenitors track each other rather closely as well. This 
is a consequence of the compensating factors that determine 
the radius of the shock during its quasi- stationary accretion 
phase. Referring to Equation 1 of Janka (12012D . these fac- 
tors (graphed for B12-WH07 and B25-WH07 in Figure^ up- 
per panel) are the PNS radius and the ^/-sphere temperature, 
whose increase increases the shock radius, and the PNS mass 
and the mass accretion rate at the shock whose increase has 
the opposite effect. Despite the variations of these factors, 
their combination during the quasi- stationary accretion phase 
(Figure [3l upper panel, solid lines) are quite similar. 

At tpb ^ 100 ms several evolutionary developments oc- 
cur in all 2D models, causing the shock trajectories of the ID 
and 2D models to diverge. First, the ratio of Tadv, the ad- 
vection timescale through the gain region (region of positive 
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Figure 3. Upper panel: Analytic shock using Equation 1 of lJankal (l20T2h 
(solid lines, dimensionless) for all simulations; and PNS radii (dotted, km), 
accretion rate at shock (dash-dotted, s"-*^), PNS mass (double-dot- 
dashed, 0.1 M0), and z/-sphere temperature (dashed, IVIeV) for the extreme 
cases, B12-WH07 and B25-WH07. Lower panel: Luminosity (solid lines, 
B s~^) and rms energies (dashed, IVIeV) for all species of neutrinos and mass 
accretion rate (dot-dashed, M© s~^) for B12-WH07. 

net neutrino heating below the shock), to Tconv, the Brunt- 
Vaisala convective growth tim escale, begins to exc eed the 
critical value Tadv/Tconv ^ 3 (iFoglizzo et alJ l2QQ6h . Con- 
vection in the gain region begins (Figure [TJ 90 ms panel), 
driven by neutrino heating and seeded by the prompt convec- 
tion episode. Neutrino-driven convection pushes the shock 
out, enlarging the gain region, and increasing Tadv Second, 
the ratio of Tadv to Theat (ratio of the thermal energy £^th to the 
net neutrino heating rate Qj,^ in the gain region) begins to ex- 
ceed unity. This signals the beginning of rising entropy in the 
gain region and an eventual shock revival ( Buras et al. 20()6a|; 
Thom pson et al.l 12005 ). Third, the SASI begins to manifest 
through low-mode distortions of the shock and oscillating pat- 
terns of the radial velocity (evident in the animated version of 
Figure[T]), causing, for example, the gain region to swell along 
the poles and contract in the equatorial regions (Figure [TJ 
150 ms panel). This distortion mode increases radv/^heat in 
the expanded regions and decreases it in the contracted re- 
gions creating outflow plumes and inflow funnels in the for- 
mer and latter regions, respectively. Thus, while the shocks of 
the ID models slowly contract after tpb ~ 100 ms, the aver- 
age shock radii of the 2D models continue to slowly increase. 
The expansion accelerates until tpb ~ 200 ms (FigureO. 
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Going fromB12-WH07 to B25-WH07, Q^^ increases - 2- 
2.5 fold, £^th increases ~ 1.5 fold, and Tadv decreases ~ 1.5 
fold, consequently radv/'?"heat ^ 1 at about the same tpb for 
all four models. On the other hand, the decrease in Tadv no- 
ticeably delays the onset of neutrino-driven convection in the 
more massive models - from tpb ^ 60 ms (B12-WH07) to 
tpb ^ 100 ms (B25-WH07). Consequently, it appears that 
neutrino-driven convection precedes the SASI in B12-WH07 
but follows it in the B25-WH07, the SASI having time to sat- 
urate in th e latter model before the onset of convection (cf. 
iMiiller et al. 2012a). B15-WH07 and B20-WH07 are inter- 
mediate cases and precedence of convection or the SASI is 
not clear. 

A Legen dre decomposition of the shock deformation (com- 
puted per iBlondin & Mezzacappal I2QQ6I) indicates that the 
growth of the SASI is dominated by the / = 2 (quadrupo- 
lar) mode in B12-WH07 and B15-WH07, while the / = 1 
(dipolar) mode is dominant in B20-WH07 and B25-WH07. 
The subsequent rise in the amplitude of shock deformation as 
the shock begins to accelerate outward is always dominated 
by the / = 1 mode. The shock deformation oscillates, with 
a period increasing with the progenitor mass, from 18 ms for 
B 12-WH07 to 30 ms for B25-WH07, though period variations 
are seen in each model. 

The neutrino luminosities Lj, and rms energies (£^r/)rms for 
all models follow a similar pattern to that of B 12-WH07 (Fig- 
ure [3j lower panel). Following the -break-out burst, the 
luminosities of all neutrino species peak between 100 and 
200 ms. The Ue-, i>e -luminosities, which arise both from the 
core and from the energy released by accreting matter, exhibit 
a more pronounced peak during the peak of the mass accretion 
rate than Lj^^^ and Lp^^, which arise more exclusively from 
the core. After 200 ms there is a rapid falloff in Lj^^ and Lp^ 
as the shock begins to accelerate outward and the mass accre- 
tion rate declines. {Ejy)rms follows the usual hierarchy, with 
energy increasing from Ue to Ue to to i^^r, the latter three 
becoming separated by only a few MeV after several hundred 
ms. The split between {Eiy^^)rms and (£^zy^^)rms is due to 
weak magnetism, which increases (decreases) the opacities of 
^/ir (^/ir)- {Eu)rms for all ncutrino species increases slowly 
as the i^-spheres contract and heat along with the PNS. 

Progenit or masses are correlated with the compactness pa- 
rameter of lQ'Connor&QttI (1201 lb (see Table [B, and the lu- 
minosities of all neutrino flavors increase with this parameter, 
nearly doubling from B12-WH07 to B25-WH07 between 50 
and 250 ms after bounce. (£^iy)rms, however, shows only a 
small increase (a few MeV) with compactness. These trends 
can be attributed to the dependence of the luminosity on 
^^-sph^z^-sph the dependence of the {Ej,)rms only on 
Tu-sph, where i^^^-sph and T^^-sph are radius and tempera- 
ture of the i^-sphere, both of which increase with the progen- 
itor compactness during the accretion shock phase (Figure [S] 
upper panel; note i^^^-gph cx i^pNs)- 

For all the models, Ljy and (^zy)rms are initially isotropic. 
By tpb = 60 ms and tpb = 100 ms for the lower and higher 
mass models, respectively, anisotropics ~10% in Li^^^p^ ap- 
pear on an / ^ 10 scale, the scale of the entropy-driven con- 
vection appearing for these models at that time. The Lj^^ 
(LpJ maxima are anti-correlated (correlated) with the con- 
vective plumes, the plumes having higher entropy but lower 
electron fraction (higher positron density). By tpb = 120 ms 
anisotropics exist at the ~ 30% level in all models on a scale 
of / = 2-3, with Lj^^ and Lp^ maxima becoming more cor- 



= sum of total energy over zones where positive 

E+gy = E"^ + overburden 

E"^Qy = E'^py + recombination energy 




time after bounce [ms] 

Figure 4. Diagnostic energy (E~^; dashed lines) versus fpb for all simu- 
lations. Dash-dotted lines (Eq^) include binding energy of overburden and 
dashed lines (£^^,rec) also include estimated energy gain from nuclear re- 
combination to ^^Ni (see definitions in text). 

related with each other and with the convective plumes and 
accretion funnels that frequently curve beneath them. By 
tpb = 200 ms, and thereafter, there is a strong / = 1 and / = 2 
component in the Li,^^p^ anisotropy, again correlated with the 
morphology of the expanding hot plumes and accretion fun- 
nels (Figure [TJ bottom four panels). Lj^^^^p^^ remain quite 
isotropic up to 200 ms after bounce, but follow, with smaller 
amplitude, the anisotropic pattern of Lj^^^p^ thereafter. 

At tpb ~ 200 ms several parameters indicate imminent 
shock revival and the onset of explosion for each model. First, 
around this time (noted as tgain in Table [B, the mass residing 
in the gain region begins to rapidly increase, which increases 
Tadv and heating efficiency ( Janka 2001; Murphv & Burrows, 
l^8l:lMarek & J anka 2009; Hanke et al. 2012).' Second, the 
^adv/^heat riscs signific antly above unity a round this time, 
signaling shock revival (iBuras et al.ll2006al: [Thompson et al.l 
2005). Third, near this time, the mean shock radius exceeds 
500 km (denoted t^oo in Table [T]). 

A useful parameter characterizing the morphology of the 
expanding shock is the shock deformation parameter, 

^ max(rshock(6^) cos 6) - min(rshock(6^) cos (9) ^ 
2 X max(rshock(6>) sin6>) 

(1) 

(IScheck et al.ll2006l) . where rshock(^) is the shock radius as a 
function of the polar angle, 6. Prolate, oblate, and spherical 
shocks have positive, negative, and vanishing values of (ighock, 
respectively. Values of (ighock for the models at tpb = 500 ms 
(Table [TJ indicate that the shock in three models is signifi- 
cantly prolate while for B20-WH07 it is nearly spherical (see 
also Figure O upper panel). 

A measure of the explosion energy frequently used is the 
diagnostic energy, (integral of total energ y, = kinetic + 
internal + gravit ational, wher e positi ve; Buras et al.l l2006al: 
ISuwa et al.ll20TQl: IMiiller et al.ll2012bh . In Figure g] we show 
E~^ for our models, as well as two other measures of the 
explosion energy: E^^ = E~^-\- overburden, and E^^ = 
E^^-\- recombination energy. The overburden is the total en- 
ergy of all negative energy zones on the grid above the inner- 
most positive energy zones plus the total energy of the off-grid 
material; the recombination energy is the energy that would 
be released if all ^He and neutron-proton pairs recombined to 
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^^Ni in the positive total energy zones. We expect that E^^ 
and E^^ bound the eventual observable explosion energy 
from below and above, respectively. 

At 500 ms all three explosion energy measures are posi- 
tive, and increasing, for B12-WH07, B15-WH07, and B20- 
WH07, but only E^ is positive for B25-WH07 due to its 
large off-grid binding energy of -0.655 B. For B12-WH07 at 
850 ms, the post-shock velocities exceed escape velocity over 
one quarter of the shock surface; the accretion rate onto the 
PNS is essentially zero; and the model is beginning to tran- 
sit from the simultaneous accretion and outflow phase to the 
wind phase. The B 12-WH07 explosion energies at 850 ms are 
E^ = 0.38 B, £;+ = 0.285 B, and £;+ ^ec = 0-32 B and ap- 
pear to be leveling off. These energies are br oadly consiste nt 
with observations for lower mass progenitors (ISmarttll2QQ9h . 

4. DISCUSSION 

We have presented results of ab initio axisymmetric core- 
collapse supernova simulations with detailed spectral neu- 
trino transport for a suite of four nonrotating models span- 
ning the mass range 12-25 Mq , through tpb = 500 ms for 
B15-WH07, B20-WH07, and B25-WH07, and tpb = 850 ms 
for B12-WH07. At 500 ms, afl models show clear indi- 
cations of developing neutrino-driven explosions, aided ini- 
tially by strong convective and SASI activity, which imparts 
pronounced prolate shock deformations in B12-WH07, B15- 
WH07, and B25-WH07, indicating that large explosion asym- 
metries are probable (in agreement with observation). 

The development of explosions in these models, taking 208- 
236 ms for the mean shock to reach 500 km, is much slower 
than in older mode ls — for exampl e, the delayed model of 
Wilson used by Woosle y et al.l (11994), for whic h t5nn ^150 
ms, or the 2D gray neutrino transport models of iHerant et al.l 
(119941) . for which tsoo ~ 75 ms. The explosion develop- 
ment times of the Chimera models reported here are com- 
parable to those of the Garching group for their models that 
use modern progenitors (tsoo ~ 250 ms for model Gil from 
iMiiller et al. ( 2012b) and tsoo ~ 230 ms for model s27 from 
iMiiller et al.l (12012a)), but significantly shorter than for their 
model s that use older proge nitors (tsoo ~ 600 ms for model 
G15 of IMiiller et al.l (l2012bl) ). This suggests a significant pro- 
genitor dependence, which makes detailed comparisons dif- 
ficult. The most pronounced difference between our models 
and those of the Garching group is the growth of the explosion 
energy. For example, once the explosion develops, E^ 'm their 
G15 model grows at a rate of 0.5 B s~^, a quarter of the rate 
in our model B15-WH07. The diagnostic energy in the Gil 
model saturates at E^ ~ 0.03, one tenth of the value of our 
model B12-WH07. The extent to which these differences are 
due to different progenitors or differences between Chimera 
and VERTEX remains to be investigated. 

Final conclusions regarding the viability of the neutrino- 
driven CCSN mechanism must wait until different groups ob- 
tain similar results from detailed ab initio simulations initiated 
from the same progenitors and carried out in 3D. Our simula- 
tions thus far support the viability of the neutrino-driven su- 
pernova mechanism for the lower progenitor masses and are 
consistent with observations. These simulations are continu- 
ing, and we will publish more complete analyses as they com- 
plete. 
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